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ABSTRACT 


A significant contributor to higher energy costs and reduced energy efficiency is the 
reactive power demand on the grid. Inductive power demand reduces power factor, 
increases energy losses during transmission, limits real power supplied to the consumer, 
and results in higher costs to the consumer. Compensating for a reactive power demand 
on the grid by providing reactive power support to the power distribution system creates 
energy efficiency gains and improves cost savings. 

One method of compensating for reactive power is by incorporating an energy 
management system (EMS) into the power distribution system. An EMS can monitor 
reactive power requirements on the grid and provide reactive power support at the point 
of common coupling (PCC) in the power distribution system in order to increase energy 
efficiency. 

The use of an EMS as a current source to achieve a unity power factor at the 
grid is demonstrated in this thesis. The power factor angle was determined using a zero¬ 
crossing detection algorithm. The appropriate amount of compensating reactive current 
was then injected into the system at the PCC and controlled using closed-loop current 
control. The process was simulated using Simulink and then validated in the laboratory 
using the actual EMS hardware. 
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EXECUTIVE SUMMARY 


A significant contributor to higher energy costs and reduced energy efficieney in 
delivering power to the consumer is the reaetive power demand on the grid. Induetive 
power demand reduces the power factor, increases energy losses during transmission, 
limits real power supplied to the eonsumer, and results in higher costs to the consumer 
due to increased power rating requirements on eleetrieal equipment. Compensating for a 
reactive power demand on the grid by providing reactive power support to the power 
distribution system creates energy effieieney gains and improves eost savings. 

One way to eompensate for reaetive power is to ineorporate an energy 
management system (EMS) into the power distribution system. An EMS can monitor 
reactive power requirements on the grid and provide reaetive power support at the point 
of common coupling (PCC) in the power distribution system in order to increase energy 
efficieney. 

The use of an EMS as a current source to achieve a unity power factor at the grid 
is demonstrated in this thesis. The power factor angle was determined using a zero¬ 
crossing detection algorithm. The appropriate amount of compensating reaetive eurrent 
was then injeeted into the system at the PCC and eontrolled using elosed-loop current 
control. The process was simulated using Simulink software and then validated using the 
actual EMS. 

A schematic of the experimental EMS’s power eleetronics is provided in Figure 1. 
Note that the EMS employs a single-phase H-bridge inverter consisting of two single-leg 
inverters. A third leg conneets a DC power supply to the inverter via a bidireetional buek- 
boost eonverter. Pulse width modulation (PWM) with unipolar voltage switching delivers 
the H-bridge gate signals that ereate an EMS output eurrent iems- A low-pass filter 
faeilitates a clean signal at the output of the inverter. Note that reactive current iems flows 
through the filter induetor. For the purpose of this thesis, the H-bridge inverter is 
eontrolled as a current source using a programmable mierocontroller that regulates the 
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polarity of the inverter output voltage and eontrols the flow of power using feedbaek 
provided by a sensor positioned at the load. 



Figure 1. EMS power eleetronics sehematie. 


The EMS output eurrent w is controlled by a reactive power compensation 
program developed in Simulink that regulates the amount of reactive current injected into 
the system. The program employs a zero-crossing detection algorithm to determine the 
power factor angle (psource at the grid. A power factor angle error correction algorithm 
eliminates any numerical error in the (psource calculation that may result from transients in 
the source current isource- Compensating Ums is generated by the EMS and controlled using 
a PI controller that adjusts the amplitude of w to drive (psource to zero. This eliminates the 
reactive power demand on the grid by achieving a unity power factor. 

The circuit used by Simulink to simulate this process is a streamlined version of 
the circuit shown in Figure 1. The simulation disregards the EMS’s power electronics by 
modeling the EMS merely as a current source. A schematic of the equivalent circuit is 
provided in Figure 2. Modeling the EMS in this fashion assumes a clean sinusoidal Ums 
signal from the H-bridge inverter without considering the PWM and closed-loop voltage 
control used to generate iems- This isolates iems in order to facilitate a focused examination 
of the Simulink-based power factor improvement methodology developed for the EMS. 
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A scenario is presented in Figure 3 in which an ohmic-inductive load is placed on 
the circuit at 0.25 sec into the simulation. This creates an overall inductive reactance in 
the circuit that demands reactive power from the grid, causing a lagging grid power factor 
of 75%. The EMS compensates by acting as a capacitive load in delivering magnetizing 
volt-amperes reactive (VAR) to the system via a capacitive iems that leads Vsource by 90°. 
The capacitive iems pulls isource in phase with Vsource- This system response is illustrated in 
Figure 3. Observe that the amplitude of iems increases over 0.30 s until iems reaches steady- 
state operation, indicating that a unity power factor is achieved at the grid. 



Time (s) 

Figure 3. Grid power factor improvement for the 0.246 H load case. 
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It is important to note that —Ums is displayed in Figures 3 to satisfy the passive 
sign convention. The passive sign convention dictates that the reference direction for 
positive current flow is into a load; however, the simulation treats the EMS as a source 
whereby w flows out of the EMS as shown in Figure 2. Hence, -Ums is presented in the 
results to show the positive flow of current into the EMS, which is consistent with the 
reference direction used to describe the positive flow of current to all other loads. 

To validate the simulation results, an experiment was conducted for the same load 
scenario. The EMS was encoded to run the Simulink program, and the ohmic-inductive 
load was placed on the EMS to replicate the aforementioned simulation scenario. The 
results of the experiment are shown in Figures 4 and 5. 

The EMS was not turned on at the start of the experiment so that the effects of the 
inductive power demand on the grid could be observed. Note in Figure 4 that the ohmic- 
inductive load creates a lagging power factor at the grid since isource lags Vsource by 
approximately 30°, which translates into an 87% source power factor. The power factor 
inconsistency between the two trials is expected since the simulation neglects the many 
electronic devices that induce reactance within the circuit.. 



Figure 4. Source voltage and current when EMS current is off. 


xviii 





































Once the EMS is turned on, capacitive power is delivered to the system via a 
leading iems, causing the phase angle of Uource to slightly lead vsource as shown in Figure 5. 
This over-compensation of reactive power to the grid occurs because the EMS does not 
currently possess a sensor on the source current and is, therefore, incapable of conducting 

closed-loop current control of 4^^. The Simulink code commands the EMS to adjust the 
amplitude of injected 4m5 ^ unity power factor is achieved, but the EMS cannot 

adjust the amplitude of 4^5 without the use of a proportional-integral (PI) controller. 

Hence, the EMS merely produces a steady 1.0 A peak (or Apk) sinusoidal 4™^ signal in 

quadrature with the source voltage, which forces the EMS to over-compensate in selling 
reactive power to the grid. Adding a source current sensor to the EMS is a planned 
hardware upgrade that should remedy the issue. 

Observe in Figure 5 the obvious presence of harmonics in the 4^5 signal. These 

harmonic signatures occur as a result of the PWM switching scheme used by the EMS to 
generate the iems signal. They are not evident in the simulation results since the circuit 
modeled in Simulink disregards the PWM associated with the EMS’s power electronics. 
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The ability of the EMS to operate as a eurrent source in compensating for a 
reactive power demand on the grid was successfully demonstrated in this thesis. The 
method utilized a power factor improvement process developed in Simulink. Specifically, 
a zero-crossing detection algorithm tracked the power factor angle between the source 
voltage and current so that an appropriate amount of reactive current can be injected into 
the system at the PCC in order to achieve a unity power factor at the grid. The process 
was simulated to predict the system’s response to capacitive and inductive power 
demands on the grid. A laboratory experiment was then conducted with the actual EMS 
to validate the process. 

The results built confidence in the ability of the EMS to compensate for a reactive 
power demand on the grid; however, improving the simulation and experiment could 
potentially facilitate a better understanding of the EMS’s capabilities. For example, the 
experiment should be repeated once a current sensor is added on the source. Also, the 
simulation could be further developed to model the EMS’s power electronics architecture 
to include the PWM and controller required to operate the H-bridge inverter. Combining 
the two aforementioned efforts into a single project might provide for a better comparison 
of the simulation and experimental results. 
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I. INTRODUCTION 


A. BACKGROUND 

The Department of Defense (DOD) has identified department-wide goals for 
increasing energy efficiency and reducing energy costs both on the battlefield and across 
its installations. In particular, the United States Marine Corps (USMC) and Department 
of the Navy (DON) seek to make 50% of all installations net-zero energy efficient and 
require that 50% of all energy consumption aboard installations come from alternative 
energy sources by 2020 [1], These goals imply that the DON and USMC must pursue 
suitable energy management practices and technologies in addition to transitioning to 
alternative fuels in order to improve energy efficiency and reduce cost at installations. 

A significant contributor to higher energy costs and reduced energy efficiency is 
the overall inductive power demand on the bulk electric grid by both DOD and civilian 
consumers. This inductive power demand reduces the power factor (PF) on the supply 
side of an electrically coupled power network, which creates numerous challenges to 
generating, transmitting, and distributing electric power to the consumer. These 
challenges are well-known and include voltage and power losses (reduced power 
efficiency) along transmission lines and in transformers, limited real (active) power 
supplied to the consumer, and higher energy and installation costs due to increased power 
rating requirements on generators and power distribution infrastructure necessary to 
support high reactive loads [2]. Transmission line length also affects cost since line losses 
vastly increase when power is transmitted over significant distances to the consumer [2]. 
This is often the case in California where power is commonly shipped from neighboring 
states. These costs are eventually passed on to the consumer. Along with ensuring energy 
security at installations [1], reducing energy costs is one of the foremost reasons why the 
DOD requires each net-zero energy installation to “produce as much energy on or near 
the installation as it consumes in its buildings and facilities” [3]. 

An example of reactive power costs charged directly to a commercial consumer in 
the UK is shown in Figure 1. Note that the installation was charged a fixed rate for a peak 
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load of 850 kVA. Real and reactive energy consumption was charged at assigned per 
unit rates; although, reactive power was charged at a much lower rate of 0.271 pence 
(0.45 cents) per kVARh (volt-ampere reactive hour). The installation’s average power 
factor during the period was relatively high at 0.949 or 94.9% [2]. 


Invoice Addrese 



SuDolv Address 


Page 3 of 4 
K 


Billing detail 


Filed Charges 


Availability charge 850 (9 £1.14/KVA 969.00 

Combe^ed HH Data Charges 20.75 

Sottlemenl agerx:y lee 1 25 

Standing charge 60.47 


TOTAL FIXED CHARGES 


Consumption Charges 

Price per 

Total 



Unit 

Unite 

Amount C 

Reactive power charge (001) 

0.00271 

613 

1 66 

Reactivo power adjustment for 31st Mar 11 

0 00271 

5789 

1568 

Unit chargeOAY RATE 

0 065654 

170482 

11192 83 

Unrt charge^GHT RATF 

0 044293 

62451 

276614 


1051 47 


TOTAL CONSUMPTION CHARGES -| 3976.31 

TOTAL ENERGY CHARGES 


Mater Reading Information 


Meter 


Unit 


Number 

Rate 

Type 

Units 

K7 

NICHT rate 

KWH 

62451 

11 

DAY RATE 

KWH 

170482 



MD 

553 



KVARH 

77481 



KWH 

232933 


E • Denotes estimated read ng, S Denotes customer reading 

Power demand and other important Information 
Power factor for this site is 0 948883 
Avord Late Payment Interest charges 


Account number 


Period of supply 

Date of invoice 

Invoice 

01 Apr 11-30 Apr 11 

03 May 11 

0027 


Figure 1. UK power bill, from [2]. 
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Commercial consumers in the United States typically are not charged directly for 
reactive power consumption. Instead, they are charged for their real power usage at 
higher peak and off-peak rates than in the UK and assessed a rebate for reactive power 
savings resulting from a high average power factor during the billing cycle. This practice 
is demonstrated by the Pacific Gas and Electric (PG&E) power bill for Naval 
Postgraduate School (NPS) in Monterey, California, shown in Figure 2. Note that NPS 
received a $458.21 rebate based on a power factor of 0.93 or 93%. The bill is split into 
two parts due to a utility rate change that occurred during the billing period. 


Details of PG&E Electric Delivery Charges 

02/18/2014 - 03/18/2014 (29 billing days) 

Service For CRNR BUTLER RttNORTH ST 
Service Agreement ID 9241321005 MAIN STATION 

02/18/2014 - 02/28/2014 

Rate Schedule E20P Service to Gusts with Max Demands of 1000 kW or 
More 

Customer Charge 11 days (g $49 28131 $542 09 

Demand Charge i 

Max Part Peak 2,251 000000 kW 3 $0.24000 204 92 

Max Demand 2,251 000000 kW 389.71000 8,290 67 

Energy Charges 

Part Peak 225.613 000000 kWh 380 09300 20.982.01 

Off Peak 215.859 000000 kWh 3 80.07734 16.694 54 

Power Factor Adjustment (3 93 00% Power Factor) -176 59 

Revenue Cycle Service Credits -1.84 

Generation Credit -28,445 36 

Power Charge Indifference Adjustnr>ent -520 03 

Franchise Fee Surcharge 278 13 

03/01/2014-03/18/2014 

Rate Schedule. E20P Service to Gusts with Max Demands of 1000 kW or 
More 

Customer Charge 18 days 3 849.28131 8887.06 

Demand Charge ^ 

Max Part Peak 2222.000000 kW 3 80.24000 331.00 

Max Demand 2222 000000 kW 3 89.71000 13.391.76 

Energy Charges 

Part Peak 295.048 000000 kWh 3 80.09451 27.884.99 

Off Peak 409,041 000000 kWh 380 07885 32.252 88 

Power Factor Adjustment (3 93.00% Rawer Factor) -281 64 

Revenue Cyde Service Credits -3.01 

Generation Credit -44.352.18 

Power Charge Indifference Adjustment -830.83 

Franchise Fee Surcharge 443.58 

Details of charges continue on next page. 

Figure 2. PG&E power bill for NPS. 
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From a comparison of Figures 1 and 2, it is apparent that energy cost savings can 
be achieved by compensating for reactive power, which directly relates to power factor 
improvement. The power factor is a measure of power efficiency and, therefore, plays an 
important role in understanding reactive power compensation. The power factor describes 
how efficiently a supply delivers real power to a load and is mathematically defined as 

PF=cos((p) ( 1 ) 

where the power factor angle (p—0^—0^ is the difference between the voltage and current 

phase angles. A unity power factor (i.e., PF= 1) implies that a sinusoidal voltage and 

current are in phase (i.e., which yields maximum real power flow since the 

relationship between power factor and real power P is given by 

P^V-I-cos(0^-0J ( 2 ) 

where voltage V and current / are the root mean square (RMS) values of the voltage and 
current. Conversely, a unity power factor implies no reactive power flow since reactive 
power Q is given by 

Q=V-I-sin(0^-0J. ( 3 ) 

Hence, achieving a unity power factor at a desired location in a circuit is the objective in 
reducing reactive power flow at the same point. 

In addition to creating energy cost savings, reactive power compensation 
increases energy efficiency gains. As previously discussed, transmission lines bleed 
reactive power more quickly than real power, which increases power losses and reduces 
power efficiency along the transmission lines—especially over long distances [4]. Also, 
most consumer loads are inductive in nature. The current of an inductance lags the 
voltage resulting in a lagging power factor, which tends to lower the system voltage [4] 
by drawing reactive energy from the system. As such, inductive loads typically cause 
voltage sags along transmission lines that limit the delivery of real power to the 
consumer. Conversely, the current of a capacitance leads the voltage resulting in a 
leading power factor, which tends to raise the system voltage [4] by delivering reactive 
energy to the system. Providing capacitive power support to the power distribution 
system, therefore, reduces power losses and improves voltage quality along transmission 
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lines [4], Transmitting reactive power also increases the amount of RMS current in the 
transmission lines, which results in additional line losses since = I^R. Hence, 

reducing reactive power demand likewise enhances the delivery of real power to the 
consumer since current otherwise wasted in power loss or transmitting reactive power is 
used to transmit real power. 

The compensation of reactive power ideally occurs as close to the inductive load 
as possible in order to eliminate the need for reactive power support from the grid by 
achieving a unity power factor [2], This typically happens at the point of co mm on 
coupling (PCC) in an electrically coupled distribution network. By providing 
compensating reactive power at the PCC, power lines from the supplier to the consumer 
are loaded only with active energy supplied by the grid [2]. Compensating to a less-than- 
unity power factor requires additional reactive power support from the grid, which further 
burdens feed-in lines to the consumer and reduces energy efficiency as a result of energy 
losses in transmission [2]. 

There are many common ways to compensate for inductive power. Most of these 
are not discussed in this thesis; however, the use of an energy management system (EMS) 
to control reactive power in a system is investigated. EMS applications in power 
electronics are becoming increasingly popular. An EMS can monitor reactive power 
requirements in an electrical circuit and inject reactive power when and where necessary 
in a circuit to increase energy efficiency within a system or limit the reactive power 
supply burden on the bulk electric grid and distributed generation (DG) sources. The 
EMS investigated in this thesis is described in more detail in the next chapter. 

B. OBJECTIVE 

The purpose of this thesis is to demonstrate the capability of an EMS to 
compensate for a reactive power demand on the grid. This is accomplished by first using 
zero-crossing detection to determine the source power factor angle (psource between the 
source voltage and current. The appropriate amount of reactive power is then injected 
into the system using the EMS as a reactive current source in order to achieve a unity 
power factor at the grid. Although the power factor improvement of a single-phase grid is 
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investigated in this thesis, the EMS can provide the same capability to three-phase 
systems. 

C. APPROACH 

In order to adequately demonstrate the ability of the EMS to facilitate reactive 
power compensation and achieve a unity power factor at the source, the first step was to 
model the EMS with zero-crossing detection capability and closed-loop current control 
using Simulink software [5], Varying loads were scheduled in the simulation to show the 
capability of the EMS to compensate for a changing reactive power demand on the grid. 
The Simulink model was streamlined in order to better isolate the zero-crossing detection 
feature of the EMS for the purpose of simplifying data collection and analysis. Certain 
model subsystems deemed unnecessary or extraneous for the purpose of this thesis were 
removed. 

Finally, a laboratory experiment was conducted using the EMS. The Simulink 
model was compiled into code that commanded the EMS to operate as a current source in 
compensating reactive power using zero-crossing detection. Experimental data was 
collected and compared to the simulation results. 

D. PREVIOUS WORK 

Much research involving the EMS discussed in this paper has been conducted to 
date. Previous Naval Postgraduate School students have added capabilities to the EMS 
over the years thereby evolving it into the current system. Recent work includes 
designing the EMS to choose between two different generators and a battery bank to 
power loads typically found at a forward operating base with the goal of improving fuel 
efficiency [6]. The ability of the EMS to store unused energy in the battery bank and 
disable generator power during periods of light loading [6] is also demonstrated. Other 
works include an analysis of the effects on microgrid power quality given different H- 
bridge inverter loading schemes [7] and using the EMS to manage peak power demand 
on a microgrid in order to improve energy efficiency and reduce energy costs [8]. 
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While the addition of a reactive power compensation capability to the EMS is 
discussed in this thesis, reactive power compensation using energy management systems 
is not a new concept. In one particular IEEE Industry Applications Magazine article, the 
authors investigate via simulation the addition of decision-making ability to a 
multifunctional single-phase voltage-source inverter (VSI) [9]. The inverter possesses the 
same smart functionalities as the EMS studied in this thesis. Specifically, the authors 
demonstrate the VSI’s ability to select from various modes of operation in response to 
certain system conditions or electricity pricing. 

In one particular case, the inverter detects a voltage sag at the grid, which is 
created by the reactive power demand of two VSI loads. The inverter then schedules 
reactive power to the PCC to alleviate the voltage sag [9]. The result is shown in Figure 
3. Note that the voltage sag is created at 0.5 s. The inverter immediately compensates, 
and the voltage sag is corrected after 1.2 s. However, the inverter continues to sell 
reactive power to the grid, causing a spike in grid voltage. The inverter’s static 
synchronous compensator function then absorbs the excess reactive power causing the 
grid voltage to stabilize [9]. 



Figure 3. Per unit grid voltage magnitude, from [9]. 


In another paper, the authors demonstrate closed-loop current control in utilizing 
DG unit “interfacing converters to actively compensate harmonics” [10]. A cost-effective 
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solution whereby fundamental and harmonic DG currents are independently controlled 
without the need for current detection in the load or voltage detection in the distribution 
system is proposed [10], The authors also use a simple closed-loop power control method 
to detect the fundamental current reference without the use of phase-locked loops thereby 
avoiding fundamental current tracking errors [10]. A similar closed-loop control 
facilitating reactive power control is later demonstrated in this thesis. 
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II. ENERGY MANAGEMENT SYSTEM 


A. FUNCTIONALITY 

The EMS utilized in this thesis is designed as an enabling technology capable of 
providing digitally-controlled “smart grid” functionality to electrical power consumers by 
facilitating power flow control to the consumer while simultaneously increasing power 
reliability and improving energy system security [11]. Depending on the situation, the 
EMS can act as a voltage source or a current source either while connected to a main AC 
source such as a bulk power grid or independent of the grid in islanded mode using 
various DG sources (photovoltaic panels, fuel cells, gas generators, batteries, etc.). 
Moreover, the EMS possesses an energy storage capability currently in the form 
of a lead-acid battery bank [11]. The following list details the major functions of 
the EMS [11]: 

• User interfacing to facilitate selective functionality. 

• Real and reactive power tracking and control. 

• Peak power management. 

• Power quality control and reliability enhancement. 

• Load scheduling and management to include non-critical load shedding 

and critical load maintenance. 

• Source power management to include DG source application. 

• Energy storage management. 

• Fault detection and correction. 

A block diagram demonstrating how the EMS might interface with its 
environment is provided in Figure 4. Note that the EMS integrates DG sources and 
storage capabilities to support various loads based on priorities established by the 
consumer. Ideally, the EMS provides power support to the grid unless a specified power 
management scenario calls for autonomous operation as previously discussed. 
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Figure 4. EMS interfacing with its environment, from [11], 


B. HARDWARE OVERVIEW 

A photograph of the EMS discussed in this paper is shown in Figure 5. Note that 
the EMS consists of power electronics framed on upper and lower printed circuit boards 
(PCB) with input/output ports and a processor cooling system. The EMS hardware is 
represented by the block diagram shown in Figure 6 and adequately explained in [11], 

The EMS’s power electronics interface with various inputs located to the 
right of the PCB is illustrated in Figure 6. The field programmable gate array (FPGA) 
development board communicates with a personal computer (PC) via an interface chip. 
The interface chip and PC couple via a joint test action group (JTAG) connector. 
Simulink code from the PC is compiled to VHDL code using Xilinx System Generator 
software [12], which is used by the FPGA to command the EMS. The FPGA 
development board is the middle board shown in Figure 5. The integrated power module 
(IPM), DC power supply, passive filters, and voltage and current sensors reside on the 
lower PCB. The IPM contains insulated gate bipolar transistor (IGBT) gate drive 
circuitry in three-phase configuration, which facilitates power flow control between the 
grid, DG sources, battery bank, and assigned loads. The voltage and current sensors 
interface with the FPGA development board via a series of analog-to-digital converters 
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located on the upper PCB. Also located on the upper PCB is a transistor-to-transistor 
interface used to command power from assigned DG sources and manage various loads. 
Finally, an LCD shows desired feedback information to the user. 



Figure 5. Photograph of the EMS analyzed in this thesis. 



Figure 6. EMS interfacing diagram. 
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C. MODELING APPROACH 

As previously discussed, the objective of this paper is to demonstrate the 
capability of the EMS to compensate for a reactive power demand on the source by 
injecting the appropriate amount of compensating reactive current into the system at the 
PCC. This capability is described in detail in the following chapter; however, it is 
important to understand the method for modeling this process in this thesis and how it 
differs from actual EMS functionality. 

A schematic demonstrating the means by which the IPM interfaces critical and 
non-critical loads with AC and DC power supplies is provided in Figure 7. The EMS uses 
a single-phase H-bridge inverter consisting of two single-leg inverters, which is preferred 
over other inverter types in high power applications [13], A third leg connects a DC 
power supply to the H-bridge inverter via a buck-boost DC-to-DC converter. A pulse 
width modulation (PWM) scheme with unipolar voltage switching delivers the H-bridge 
gate signals [7] and reduces the effects of harmonics in the output voltage at the 
switching frequency [13]. A low-pass EC filter (EPF) facilitates a clean AC signal at the 
output of the inverter. The EPF is shown in Figure 7. Note that the filter capacitor is 
marked by Vac, and reactive current w flows through the filter inductor. 



Figure 7. EMS power electronics circuit schematic. 
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Recall that the H-bridge inverter can be controlled as a current source or a voltage 
source. For the purpose of this thesis, the H-bridge inverter is controlled as a current 
source using a programmable microcontroller. The microcontroller provides a number of 
control features to the EMS. In particular, it regulates the polarity of the output voltage 
from the H-bridge inverter by controlling the H-bridge inverter’s IGBT switches and 
controls the flow of power using feedback provided by a sensor positioned at the load. 
This control scheme is shown in Figure 7. 

As a current source, the EMS produces an output current w, which is used to 
inject compensating reactive power at the PCC. A particular method for controlling igms in 
providing reactive power support to the AC grid is demonstrated in this thesis. This 
method involves using a Simulink-based zero-crossing detection algorithm to determine 
the power factor angle (psource at the grid and, subsequently, changing the amplitude of the 
compensating reactive current w to bring (psource to zero. Achieving a unity power factor 
at the source eliminates reactive power demand on the grid. 

To effectively model this process, an equivalent of the circuit shown in Figure 7 is 
simulated using Simulink software. The reason for using an equivalent circuit is to 
demonstrate the EMS as a constant current source in order to examine the use of Ums in 
compensating for a reactive power demand on the AC grid. A schematic of the idealized 
circuit is illustrated in Figure 8. Modeling the EMS in this fashion assumes a clean 
sinusoidal Ums signal from the H-bridge inverter without having to consider pulse width 
modulation and associated closed-loop voltage control to produce the w necessary to 
achieve a unity power factor at the grid. 
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Notice from Figure 8 that the EMS block in Figure 7 is effectively eliminated 
from the idealized circuit. All circuitry to the left of the 12 gF filter capacitor in Figure 7 
is not considered in the schematic in Figure 8 except for the flow of w from the EMS to 
the PCC. This idealized modeling scheme facilitates a more focused examination of the 
power factor improvement methodology developed for the EMS using Simulink, which is 
explained in detail in Chapter 111. 

The experimental EMS of course remains as shown in Figure 7 and, therefore, 
employs pulse width modulation to create the H-bridge gate signals that generate a 
voltage to control The power factor improvement methodology developed in 
Simulink still applies to the experiment, but no proportional-integral (PI) controller was 
implemented in the lab to correct Ums because the EMS hardware currently does not have 
a sensor on the source current. Adding a sensor to the source current is a planned upgrade 
to the system. Due to this restraint, the experiment was conducted without the use of 
closed-loop current control, but the EMS still demonstrated the capability to generate 
compensating reactive current in response to a reactive power demand on the grid. The 
results are explained in further detail in Chapter IV. 


14 



III. COMPUTER SIMULATION 


A. OVERVIEW 

A computer simulation of the idealized circuit shown in Figure 8 was created 
using Simulink. The simulation design is based on the circuit shown in Figure 9, which is 
the same circuit presented in Figure 8 except with two defined variable loads tested 
both in simulation and experiment. Load 1 is a purely resistive 85.7 Q. load. Load 2 is a 
0.246 H inductive load with an estimated internal resistance of 5 O. The purpose of the 
two loads is to demonstrate the capability of the EMS to compensate for a varying 
reactive power demand on the grid. 



The circuit is mathematically simulated by the circuit component and subsystem 
diagram in Figure 10. The 12 pF filter capacitor and 400 pH filter inductor with 10 mQ 
resistance are clearly displayed in the diagram as independent of the subsystems. A 
constant 110 Vrms, 60 Hz sinusoidal input signal from the grid precedes the filter inductor. 
The Load Switching subsystem contains both variable loads, which are controlled by an 
electronic switch that adds inductance to the system at a predetermined time in the 
simulation. The Power Factor Correction subsystem contains the zero-crossing detection 
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and power factor angle error correction algorithms used to continuously track phase 
angles of vsource and isource (and 6 , respectively) and subsequently determine and 

^source ^source 

correct power factor angle (psource in order to achieve a unity power factor at the grid. The 
Power Factor Correction subsystem also employs closed-loop current control, which 
drives (psource to zero by changing the amplitude of iems as it flows to the PCC. 
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Figure 10. Simulink circuit component and subsystem diagram. 


To keep the model as visually clean as possible, all circuit data regarding output 
variables 1 through 11 (depicted in Figure 10) is sent to the top-level block diagram 
shown in Figure 11. There it is collected along with simulation time and sent to the 
MATLAB “workspace” for graphical analysis. The MATLAB file scripts containing the 
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system’s initial conditions, discrete component values, and code used to create the plots 
in this thesis are presented in the Appendix. 
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Figure 11. Simulink top-level block diagram. 


B. LOAD SWITCHING 

As previously stated, load switching is used to demonstrate the capability of the 
EMS to compensate for a varying reactive power demand on the grid. The resistive and 
inductive loads shown in Figure 9 are simulated using the load switching subsystem 
presented in Figure 12. At the start of the simulation, only the 85.7 Q. load is active, 
which results in an overall capacitive reactance in the circuit due to the dominating 
presence of the 12 pF low-pass filter capacitor in the circuit. The result is an overall 
capacitive power demand on the source that is quickly compensated by the injection of 
inductive iems at the PCC. 

At 0.25 s into the simulation, the 0.246 H inductive load with 5 Q internal 
resistance is added to the circuit. The load switching operation is controlled by the load 
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step and Switch 1 function blocks shown in Figure 12. This new load scheme creates an 
overall inductive power demand on the grid that models the power demand of a typieal 
reactive power consumer. A capacitive w is accordingly injected at the PCC to correct 
the lagging power factor produced at the grid by the inductive load. The results are 
presented in Section D of this chapter. 


Load 1 R (85.7 Ohms) 


v_cap 

Cl> 



Load 2 R (5 Ohms) 

Figure 12. Simulink load switching subsystem diagram. 


C. POWER FACTOR CORRECTION 

The method used to create a unity power factor at the source is conceptually 
simple. The power factor correction flow chart shown in Figure 13 illustrates the process. 
First, the source power factor angle (psomce is calculated using zero-crossing detection. As 
(Psource is gradually driven to zero via closed-loop control, it is passed through a simple 
first-order, low-pass filter transfer function. The LPF filter softens the change in (psource 
over time, which assists the PI controller in smoothly adjusting reference current 

amplitude . The EMS concurrently generates a continuous 1.0 Apt sinusoidal current 
in quadrature with the sinusoidal source voltage that is multiplied with to create the 
sinusoidal EMS reference current 4^^. The referenee eurrent is designed to lead 


18 

























Vsource by 90° whcn isource lags Vsource SO that 4m* impTOves the power factor at the grid 

T* 

when injected into the system. Alternatively, the reference current amplitude will be 

negative so that the injected lags vsource by 90° in order to improve a leading grid 
power factor. 

Conceptually, 4m* is then compared to the measured EMS output current 4m* 
the flow chart. The difference between the two waveforms defines the error that is 
corrected by another PI controller, which directs the inverter to generate a new 4m* '^i^ 

changes in the PWM scheme until 4m* ^ Cm* ■ ibis point, the system makes no further 

changes to 4m* since a unity power factor is realized at the grid. Recall that the 

experimental EMS currently lacks a sensor at the source, and it is not yet possible to use 
closed-loop control to makes changes to the actual PWM in the laboratory. The results of 
the experiment are discussed in Chapter IV. 



Figure 13. Power factor correction flow chart. 


Recall from Figures 8 and 9 that the Simulink model does not account for PWM 
in the EMS and, therefore, does not create an 4m* signal. The simulation instead 

simplifies the power factor correction process by assuming that 4m* ^ Cm* thus, only 

adjusts as necessary to correct the source power factor. Accordingly, feeds 

continuously into the PCC and is regulated via the first PI controller shown in Figure 13 
as (psource chauges over time. 
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The function of the power factor correction diagram shown in Figure 14 is to 
implement the power factor correction flow chart in Figure 13 using Simulink. The 
components and subsystems that comprise the power factor correction block are 
illustrated by the Simulink power factor correction subsystem diagram in Figure 14. Note 

that the diagram replicates all flow chart functions leading to the creation of the 

signal. The PI controller and quadrature reference current signal are contained within the 
closed-loop current control subsystem. All parameter data is sent to the output of the 

power factor correction diagram to include 4^^, which is then injected back into the PCC 
as demonstrated in Figure 10. 
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Figure 14. Simulink power factor correction subsystem diagram. 


1. Zero-Crossing Detection 

In order to determine the source power factor angle (psource, the EMS uses zero¬ 
crossing detection algorithms to determine phase angles 0^ and 6: . These 

^source ^source 

algorithms are located in the source voltage and source current Zero-Crossing Detection 
subsystems illustrated in Figure 14. The algorithm for determining 6^ is less complex 
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since Vsource remains constant. The Simulink Vsource zero-crossing detection algorithm for 


determining 6 is shown in Figure 15. 

^source 


^^ double 

Q3 

v_source 



Rise (Voltage) 


Figure 15. Simulink source voltage zero-crossing detection diagram. 


The algorithm detects the positive rise of Vsource by determining if the present 
sinusoidal voltage input value is positive and the previous value is non-positive. If a 
positive rise is detected, the Detect Positive Rise function block outputs a Boolean 
TRUE. Otherwise, the block outputs a Boolean FALSE. A TRUE statement signifies that 
the signal crossed the time axis while rising, which occurs once every cycle for a 
sinusoidal signal. The Integrator function block integrates over time t of the source 
voltage signal frequency from zero to a reset time of period T and then repeats. The 
integrator is programmed to reset when the Detect Positive Rise function block outputs a 
Boolean TRUE (i.e., detects a positive rise across the time axis). The result gives the 
phase angle by the relationship 

^ £ Itt fdt (4) 

where the period T is the time between successive positive rises across the time axis and 
/is the signal frequency, which is typically 60 Hz for the grid voltage. 

The EMS computes 0- using an algorithm similar to that of 6 , except the 

''source ^ ^ ^source ^ 

isource zero-crossing detection algorithm compares each integrator output value of to 
the quantity 1. 671 . The algorithm ensures that the present value of ^ is greater than 

1 .671 during a positive rise before reporting a Boolean TRUE to the Integrator function 
block. This prevents the integrator from prematurely resetting in the event that signal 
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transients or disturbances in Uource cause the signal to unexpectedly perform a zero¬ 
crossing. This process is demonstrated by the Simulink diagram in Figure 16. 
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Figure 16. Simulink source current zero-crossing detection diagram. 




Figure 17. Source voltage and current phase angle plots for both load cases. 


Plots showing the change in 0^ and 0: over time for both the resistive and 

^ ^ ^source ^source 

inductive load cases are presented in Figure 17. The graphic only shows the first 0.35 s of 
the simulation for the purpose of visual acuity. As previously mentioned, an overall 
capacitive reactance initially exists in the circuit before adding the inductive impedance 

to the load resulting in a leading power factor since isource leads Vsource- Once the inductive 
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load is introduced, the circuit experiences an overall inductive reactance, which causes 
isource to lag Vsource and results in a lagging power factor. Notice that the simulation 
gradually brings isource in phase with Vsource for both load cases since the simulation is 
running the power factor correction algorithm in the example. 

2. Power Factor Angle Error Correction 

The phase difference between 0^ and 6 is then computed to determine the 

^source ^source 

source power factor angle (psource by subtracting the two quantities via the Add5 block 
shown in the Power Factor Correction subsystem in Figure 14; however, computing 
(Psource in this fashion results in numerical errors. The errors are caused by the 
asynchronous resetting of the Integrator function blocks in Figures 15 and 16 in 

computing 0^ and 6 , respectively, using zero-crossing detection. 

This mathematical glitch is best understood from Figure 17. In the inductive load 
case, for example, isource lags Vsource allowing the 0y^ integrator to wrap from 2?! to zero 

before the ^ integrator. The result is a momentarily large cpsource that abruptly declines 

as (psource approaches zero. These errors are displayed as spikes in the power factor 
calculation plot shown in Figure 18 and must be eliminated before closed-loop current 
control can operate to pull isource in phase with Vsource- 



Time (s) 

Figure 18. Power factor calculation with numerical integration error. 
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The Power Factor Angle Error Correction subsystem shown in Figure 19 
eliminates the numerical error in (psource simply by adding In to (psource if (Psource is less than 
n or subtracting 2n from cpsource if (Psource is greater than n. This arithmetical process is 
possible since the calculated error in ^source cannot be greater than ±2;! simply because 
(Psource caunot mathematically be more than ±2;! as demonstrated by the error spikes in 
Figure 18; hence, numerical error nears ±2;! as (psource approaches zero. Note in Figure 18 
that the power factor angle is plotted in degrees instead of radians. 



Figure 19. Simulink power factor angle error correction diagram. 


The corrected power factor angle is then passed through a simple first-order low- 
pass filter transfer function as seen in Figure 14. The LPF softens the change in (psource 
over time, which assists the Closed-Loop Current Control subsystem in smoothly 

adjusting the reference current amplitude . Plots of the corrected and filtered power 

factor angle are shown in Figure 20. Note from the plots that cpsource < 0 when the circuit 
is capacitive (0 s < time < 0.25 s) and cpsource > 0 when the circuit is inductive (time > 0.25 
s). The load change is scheduled at precisely 0.25 s; however, the zero-crossing detection 
algorithm does not detect a change in 9^ until the end of the next zero-crossing event, 

^ ‘‘source ^ 

which occurs one full period (1/60 s) later. 
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Figure 20. Corrected and fdtered power factor angle plots. 


3. Closed-Loop Current Control 

The function of the Simulink Closed-Loop Current Control subsystem is to 
modify EMS output current by changing its amplitude . The closed-loop current 
control block diagram is shown in Figure 21. The diagram mathematically replicates the 
generation process described by the flow chart in Figure 13. 
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Figure 21. Simulink closed-loop current control diagram. 
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T* 

EMS current amplitude is produced as the output of the PI eontroller. The 
controller uses negative feedback to drive (p to zero. The fed-baek output signal is 
reactive current 4ms» which is created by multiplying with a 1.0 Apk eontinuous 

sinusoidal waveform generated in quadrature with Vsource- The fed-back output signal thus 
has the form 

Cs=Csin(^,_-^) (5) 

T* 

where amplitude is adjusted by the PI controller to bring isource in phase with Vsource- 

.* 

Changing the phase angle of isource by adjusting the amplitude of is possible 
by Kirehhoff s Current Law (KCL). This is apparent when observing the flow of eurrent 
through the PCC of the eireuit in Figure 9. Given that 4ms ^ hms the simulation 
cireuit, the KCL equation for isource is 

i =i, ,+i —i (6) 

source load cap ems \ / 

where positive 4ms deseribes the flow of EMS eurrent into the PCC, henee the minus 
sign in Equation (6). 

Plots of and (psource versus time are displayed in Figure 22. The plots 
demonstrate the relationship between and (psource —speeifically how EMS current 
amplitude 1^^ changes aceordingly to pull isource in phase with Vsource (i-c., foree (psource to 
zero) thereby aehieving a unity power factor at the grid. Also note that approaches 
unity as (psource approaches zero. This occurs because no further ehange to is required 

when (Psource = 0. Observe that can be negative or positive depending on whether 
isource Icads or lags Vsource, respeetively. 
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Figure 22. Changing EMS current amplitude to bring the source current in 
phase with the source voltage in order to achieve a unity power factor. 


D. RESULTS 

The simulation results are graphically displayed in Figure 23 and Figure 24. The 
figures demonstrate the system responses in terms of Vsource, hource, and iems for the 
resistive and inductive load cases exhibited in Figure 9. The results shown in Figure 23 
are in response to the 85.7 Q. resistive load placed on the EMS at the start of the 
simulation. At 0.25 s, the 0.246 H inductive load is added to the circuit. The system 
response to the inductive load is illustrated in Figure 24. In both figures, Vsource is scaled 
down by a factor of 50 in order to more adequately contrast the three signals. The voltage 
Vsource Still rcprcsents the 110 Vrms voltage supply provided by the grid. 

At the start of the simulation, the circuit impedance is dominated by the 12 pF AC 
filter capacitor despite the active 85.7 Q load. Hence, the overall circuit reactance 
is capacitive resulting in a leading grid power factor where initially (psource = “21° as 
illustrated in Figure 23. The leading power factor implies that the filter capacitor is 
selling reactive power to the grid. The EMS compensates by generating an inductive iems 
that lags Vsource by 90°. This EMS effectively absorbs reactive power from the system in 
order to bring the hource signal in phase with Vsource- Notice that the amplitude of iems 
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gradually increases to 0.71 Apk over 0.24 s until Ums reaches steady-state conditions, 
which indicates that a unity power factor is achieved at the grid. 
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Figure 23. Grid power factor improvement for the purely resistive load case. 


At 0.25 s, the 0.246 H inductive load is added to the system generating an overall 
inductive reactance in the circuit. The inductance demands reactive power from the grid 
causing a lagging grid power factor. The results in Figure 24 show isource initially lagging 
Vsource by 41° upon activation of the inductive load, which results in a grid power factor of 
just over 75%. The EMS compensates by acting as a capacitive load in delivering 
magnetizing VARs to the system. The EMS injects capacitive w at the PCC where iems 
leads Vsource by 90°. By delivering reactive power to the system, the EMS gradually brings 
isource into phase with vsource- In this case, the amplitude of Ums increases to about 1.01 Apk 
over 0.30 s until iems reaches steady-state conditions as a unity power factor is achieved at 
the grid. At this point, the EMS makes no further adjustment to Ums- 
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Figure 24. Grid power factor improvement for the inductive load case. 

It is important to note that - w is used to represent EMS current in the simulation 
results shown in Figure 23 and Figure 24 in order to satisfy passive sign convention. 
Passive sign convention dictates that the reference direction for positive current flow is 
into a passive circuit component (i.e., a load). The load current iioad and filter capacitor 
current icap shown in Figures 8 and 9 are properly referenced using passive sign 
convention. This is also evidenced by Equation (6) using KCL; however, the simulation 
treats the EMS as an active circuit component (i.e., a source) whereby Ums flows from the 
EMS to the PCC as shown in Figure 9. Hence, -iems is given in Figures 23 and 24 since it 
portrays the positive flow of current into the EMS, which is consistent with the reference 
direction used in the simulation to describe positive current flow through all other loads. 
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IV. LABORATORY EXPERIMENT 


A. SETUP 

The experiment was set up in the laboratory using the EMS hardware described in 
Chapter 11. Recall that the experimental EMS design incorporates the power electronics 
circuitry shown in Figure 25, which involves the use of PWM to generate the IGBT gate 
signals of the H-bridge inverter. The main electric grid provided the 60 Hz, 110 Vrms AC 
source input to the EMS. Clamp-on probes were used to take source voltage and current 
measurements, which were captured using an oscilloscope. A 72 Vdc battery bank 
powered the 200 Vdc bus via the DC-DC boost converter shown in Figure 25. The DC 
bus provided the DC input voltage to the H-bridge inverter. 

All discrete component values for the circuit in Figure 25 are given in Table 1. 
Also, a digital photograph of the EMS under test is shown in Figure 26. 



Figure 25. Experimental EMS power electronics circuit schematic. 
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Table 1. 


Discrete component values for the circuit in Figure 25. 


Device 

Value 

AC grid 

110 Vnns, 60 Hz 

Vbatt 

72 Vdc 

Lin 

200 pH 

Ldc 

300 pH 

Lems 

200 pH 

Cin 

12 pF 

Cdc 

500 pF 

Cbus 

2niF 

Load 1 

85.7 Q 

Load 2 

0.246 H 



Figure 26. The EMS under test in the laboratory. 
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The oscilloscope, PC with Simulink and Xilinx software, and variable circuit 
loads are not shown in Figure 26; however, the variable load panels are shown in the 
photographs in Figure 27. The circuit loads were controlled by two load panels that 
connected to the EMS. Each panel consists of three parallel sets of three parallel resistors 
(300 Q, 600 Q, and 1200 Q.) or inductors (0.8 H, 1.6 H, and 3.2 H) to support up to three- 
phase load applications. The variable load boxes were accordingly configured to facilitate 
the experiment. This setup can be seen in Figure 27 where the appropriate resistor and 
inductor switches were activated to provide the 85.7 Q. and 0.246 H parallel loads used in 
the experiment. 



Figure 27. Variable load panels used in the experiment. 
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B. PROCEDURE 


Implementation of the experiment was a straightforward process. Upon 
completing the lab setup, the Simulink simulation was compiled into VHDL code on a 
PC and uploaded to the FPGA development board on the EMS using Xilinx software. 
The VHDL code enabled the EMS to execute the power factor correction methodology 
developed in Simulink. 

The purpose of the experiment was to validate simulation results for the 
compensation of an inductive power demand on the grid. Thus, both the 85.7 Q resistive 
load and 0.246 H inductive load were effective at the start of the experiment in order to 
emulate circuit conditions similar to those modeled with Simulink upon activation of the 
inductive load at 0.25 s into the simulation. 

EMS functionality permitted the manual activation and deactivation of ; 

therefore, was not immediately turned on so that the effects of the inductive power 
demand on the grid could be observed and recorded using the oscilloscope. Once initial 
measurements were collected, was activated, and the resultant waveforms were 
captured. The resultant experimental data is presented in Figures 28-31. 

Since the EMS did not possess a current sensor at the source as previously 
mentioned, the microcontroller was incapable of measuring 6. and, therefore, could 

^source 

not conduct closed-loop control of reference current 4^^. The EMS could only generate a 

single 4ms signal from the H-bridge inverter. Consequently, the EMS was merely able to 
produce the 1.0 Apk continuous quadrature current waveform shown in Figure 13 to 
compensate for a reactive power demand on the grid. Adjustments to based on 
changes to 9. were not possible with the current hardware configuration. However, 

^ ''source ^ 

the purpose of the experiment was to demonstrate the potential of the EMS to compensate 
for a reactive power demand on the source using the method tested in simulation. This 
capability is confirmed by the experimental results. 
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C. RESULTS 

The ohmic-inductive circuit load initally creates an overall inductive reactance in 
the circuit that results in a lagging power factor at the grid as expected. This is shown by 
the source voltage and current waveforms in Figure 28 where isource lags vsource by 
approximately 30°, which translates into an 87% source power factor. This power factor 
is slightly improved over the 75% power factor initially realized in simulation; however, 
this discrepancy is expected since the simulation does not account for the many electronic 
devices that induce reactance throughout the EMS circuitry. 

Notice also in Figure 28 that no reactive power compensation occurs since 4^5 is 

not yet activated. There also appears to be small negative DC voltage component in the 
source voltage of about -0.34 V. This DC offset is a consequence of calibration error in 
the EMS voltage sensor since the equipment was not properly calibrated before the 
experiment was conducted. Also of note, the voltage and current waveforms shown in 
Figure 28 were created in MATLAB using data imported from the oscilloscope. The 
oscilloscope display of the very same waveforms is provided for reference in Figure 29. 



Figure 28. Source voltage and current when EMS current is off (MATLAB). 
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Figure 29. Source voltage and current when EMS current is off (oscilloscope). 


Once is activated, capacitive power is delivered to the system causing the 

phase angle of Uource to shift until isource slightly leads Vsource- This over-compensation of 
reactive power to the grid occurs because the EMS is currently incapable of employing a 
proper closed-loop current control scheme. The EMS merely produces a constant 1.0 Apk 

sinusoidal signal in quadrature with the source voltage. The amplitude of 4^^ does 

not change without a PI controller, and a leading power factor occurs at the grid. This 
situation is illustrated in Figure 30. A proper closed-loop current control method would 

otherwise adjust the amplitude of 4^^ so as to pull isource in phase with Vsource- 

Notice in Figures 30 and 31 that is again used to represent EMS current. 

This is because the EMS executes the same power factor correction algorithm used by the 
simulation whereby the EMS is treated as a current source in accordance with passive 
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sign convention; thus, presenting in the experiment results serves to characterize 

the EMS as a circuit load, which ensures consistency in referencing the direction of 
positive flow of current from the PCC to the loads. 

The presence of harmonics in the 4^^ signal is also obvious from the 

waveforms shown in Figures 30 and 31. The harmonics stem from the PWM switching 
scheme used by the EMS to generate iems- These harmonic signatures are not evident in 
the simulation results since the circuit modeled in Simulink disregards PWM associated 
with the EMS’s power electronics. 



0.015 


0.02 


Figure 30. Source voltage and current when EMS current is on (MATLAB). 
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Figure 31. Source voltage and current when EMS current is on (oscilloscope). 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Reactive power compensation and power factor improvement are not new 
engineering concepts. Controlling the generation, transmission, and distribution of 
reactive energy in delivering quality power to the consumer is essential to increasing 
energy efficiency and reducing energy costs. As a result, various methods for achieving a 
unity power factor at the source have been developed and improved over time. 

A particular means of compensating for a reactive power demand on the grid was 
examined in this thesis whereby an additional capability to a particular EMS was 
proposed that enabled the EMS to operate as current source in compensating for a 
reactive power demand on the grid. This was accomplished by first using a Simulink- 
based zero-crossing detection algorithm to determine the power factor angle between the 
source voltage and current. The appropriate amount of reactive current was subsequently 
injected into the system and adjusted using a closed-loop current control scheme that 
brought the source current in phase with the source voltage thereby eliminating any 
reactive power demand on the grid. A Simulink model of the process was initially 
developed in order to forecast the system’s response to both capacitive and inductive 
power demands on the grid. The process was then confirmed in a laboratory using the 
actual EMS. 

It is important to remember that the Simulink model was simplified so as to 
isolate the reactive power compensation process for analysis by specifically neglecting to 
model the EMS’s power electronics systems. The PWM scheme and related PI control of 
the H-bridge IGBT gate signals in particular were disregarded. The EMS was thus 
modeled as a constant current source. This eased the complexity of the simulation design 
but created disparities between model functionality and actual EMS operation, which 
caused some dissimilarity between simulation and experimental results. 

For example, no harmonics were observed in the simulation representation of 4^5 

since the effects of the EMS’s power electronics were disregarded by the circuit modeled 
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in Simulink. Likewise, additional circuit reactance associated with discrete electronic 
components was not observed in the simulation and, therefore, did not affect the reactive 
power demand on the grid. Nevertheless, the results verified the ability of the EMS to 
apply the Simulink-based power factor correction algorithm in compensating for a 
reactive power demand on the grid. 

Moreover, the EMS hardware did not possess a current sensor at the source, so 
no closed-loop control of was implemented in the laboratory. The EMS was 

consequently unable to adjust the amplitude of 4^^ in correcting the source power factor 

angle. The EMS still demonstrated the ability to improve the power factor at the source, 
but a unity power factor could not be achieved without the use of closed-loop current 
control. An EMS hardware upgrade is currently planned that will improve the 
functionality of the EMS. 

B. RECOMMENDATIONS 

The results presented in this thesis built confidence in the ability of the EMS to 
compensate for a reactive power demand on the grid; however, improving the simulation 
and experiment could potentially facilitate a better understanding of the EMS’s 
capabilities. It would be advantageous, for example, to repeat the experiment once a 
source current sensor is installed in the EMS in order to fully validate the reactive power 
compensation process designed in Simulink. Also, the simulation could be further 
developed to simulate the EMS’s power electronics architecture, which includes the 
PWM and controller required to operate the H-bridge inverter. Combining the two 
aforementioned efforts into a single project might enable a better comparison of the 
simulation and experiment results. 

Additionally, future research into other reactive power compensation methods 
involving the EMS would add depth to its functionality and facilitate a comparative 
analysis of the different methods. Investigating the system’s response to various reactive 
power demand scenarios might also assist in identifying the advantages and limitations of 
each method. These results could then be used to develop a decision-making algorithm 
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that enables the EMS to choose an optimal power factor improvement process when 
faced with varying power supply and demand situations. 

One example of a potential alternative to the reactive power compensation 
method presented in this thesis is the use of a phase-locked loop (PEL) control scheme to 
detect and match the source current and voltage phase angles. While developing a PEL 
control scheme may be a more complex design effort, it could potentially alleviate the 
numerical integration error associated with zero-crossing detection when signal transients 
are present in the source current. 
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APPENDIX. MATLAB M-FILES 


A. SIMULATION INITIAL CONDITIONS FILE 

% Michael Prato 

% Reactive Power Compensation Using an EMS 
% Initial Conditions File 


clear all; 



clc; 



close all; 



V_source=120’^sqrt (2) ; 

o, 

o 

Source voltage (V peak) 

f req=2 ’^pi’^ 60 ; 

o, 

o 

Source voltage frequency (rad) 

Kp=2/4; 

o, 

o 

Proportional gain of the PI controller 

Ki=100/4; 

o. 

Integral gain of the PI controller 

R_load=85.714; 

o. 

Resistive load (ohms) 

L_load=0.246; 

o, 

o 

Inductive load (H) 

R_o f_L_1o ad= 5; 

o, 

o 

Resistance of the inductor (ohms) 

Rin=0.01; 

o, 

o 

Internal resistance (ohms) 

L_fil=400e-6; 

o. 

Filter inductor (H) 

C_fil=12e-6; 

o. 

Filter capacitor (F) 

tstep = 0.5e-5; 

o, 

o 

Step size (sec) 

tstop = 0.6; 

o, 

o 

Sim length (sec) 


B. SIMULATION OUTPUT PLOT FILE 


% Michael Prato 

% Reactive Power Compensation Using an EMS 
% Output Plot File 

% Replicates oscilloscope waveforms based on output data files 
collected from the oscilloscope 

v_cap=data_out(:,1); 
i_source=data_out(:,2); 
i_ems=data_out(:,3); 
i_load=data_out(: , 4) ; 
filtered_PF_angle=data_out(: , 5) ; 
v_source=data_out(: , 6) ; 
i_ems_amplitude=data_out {:,!)} 
i_phase_angle=data_out(:^ 8); 
v_phase_angle=data_out(:,9); 
calculated_PF_angle=data_out(:,10); 
corrected_PF_angle=data_out(:,11); 

figure( 'name'Reactive Power Compensation'); 
hold on; 

plot(time^ v_source/50^ 'b' , 'linewidth' ,1.5) ; 
plot(time,i_source, 'r' , 'linewidth' ,1.5) 
plot(time,-i_ems, 'g' , 'linewidth' ,1.5) 

xlabel( 'Time (s) ' , ' FontSize' ,20) ; ylabel( 'Amplitude' , 'FontSize' ,20); 
xlim([0 tstop]); 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '- i_e_m_s (A)' ); 
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grid on; 
hold off; 

figure( 'name' ^ 'Grid Power Factor'); 

plot(time^ cos(calculated_PF_angle), 'linewidth' , 2 )} 

xlabel( 'Time (s)' ^ 'FontSize' ,20); ylabel( 'Grid Power 

Factor' ' Font Size' 20) ; 

xlim([0 tstop]); ylim([0.75 1.02]); 

grid on; 

figure( 'name'Reactive Power Compensation by Achieving a Unity PF, 

Ohmic Load' ); 

subplot (2^1,1); grid on; 

hold on; 

plot(time,v_source/50, 'b' , 'linewidth' ,2); 
plot(time,i_source, 'r' , 'linewidth' ,2) 
plot(time,-i_ems, 'g' , 'linewidth' ,2) 

%title('Purely Resistive Load (85.7 \Omega)','FontSize',20); 
xlim([0 0.25]); 

xlabel( 'Time (s)' , 'FontSize' ,20); ylabel( 'Amplitude' , 'FontSize' ,20); 
legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '- i_e_m_s 

(A)' , 'location' , 'south' ) ; 
subplot(2,1,2); 

plot(time,cos(calculated_PF_angle) ,' linewidth' , 2) ; grid on; 
xlabel( 'Time (s)' , 'FontSize' , 20); ylabel( 'Grid Power 
Factor' , 'FontSize' ,20); 

%title( 'Power Factor','FontSize',20) ; 
xlim([0 0.25]); ylim([0.93 1.005]); 

figure( 'name' ,' Reactive Power Compensation by Achieving a Unity PF, 
Inductive Load' ); 
subplot(2,1,1); grid on; 
hold on; 

plot (time,v_source/50, 'b' , 'linewidth' , 2) ; 
plot(time,i_source, 'r' , 'linewidth' , 2) 
plot(time,-i_ems, 'g' , 'linewidth' , 2) 

%title('Inductive Load (0.246 H) Added to the Circuit at t = 0.25 
s','FontSize',20); 

xlabel( 'Time (s)' , 'FontSize' ,20); ylabel( 'Amplitude' , 'FontSize' ,20); 
xlim([0.25 tstop]); 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '- i_e_m_s 

(A)' , 'location' , 'south' ) ; 
subplot(2,1,2); 

plot(time,cos(calculated_PF_angle) ,' linewidth' , 2) ; grid on; 
xlabel( 'Time (s) ',' FontSize' , 20); ylabel( 'Grid Power 
Factor' , 'FontSize' ,20); 

%title( 'Power Factor','FontSize',20) ; 
xlim([0.25 tstop]); ylim([0.75 1.02]); 

figure( 'name' , 'Phase Angles' ) 

plot (time, v_phase_angle’^180/pi, ' b' , ' linewidth' ,1.5) ; 
hold on; 

plot (time, i_phase_angle’^180/pi, ' g' , ' linewidth' ,1.5) ; 
hold off; 
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xlabel( 'Time (s) ' ^ ' FontSize' ^16) ; ylabel( 'Phase Angle, \theta 
(deg)' , 'FontSize' , 16) ; 

legend( '\theta_v_{_s_o_u_r_c_e}' , '\theta_i_{_s_o_u_r_c_e}' , 'FontSize' ,1 
6) ; 

xlim([0 tstop]); ylim([0 410]); 
grid on; 

figure( 'name' ,' Phase Angles from t=[0, 0.6]sec' ) 
subplot(2,1,1); grid; 
hold on; 

plot (time, v_phase_angle’^180/pi, ' b' , ' linewidth' ,1.5) ; 
plot (time, i_phase_angle’^180/pi, ' g' , ' linewidth' ,1.5) ; 
hold off; 

title( 'Resistive Load (85.7 \Omega)' , 'FontSize' ,20); 
xlim([0 0.25]); 

xlabel( 'Time (s)' , 'FontSize' ,16); ylabel( 'Phase Angle, \theta 
(deg)' , 'FontSize' , 20) ; 

legend ( ' \theta_v_{_s_o_u_r_c_e }' , ' \theta_i_{_s_o_u_r_c_e }' ) ; 
subplot(2,1,2); grid; 
hold on; 

plot (time, v_phase_angle’^180/pi, ' b' , ' linewidth' ,1.5) ; 
plot (time, i_phase_angle’^180/pi, ' g' , ' linewidth' ,1.5) ; 
hold off; 

title (' Inductive Load (0.246 H) Added to the Circuit at t = 0.25 
s' , 'FontSize' ,20); 
xlim([0.25 tstop]); 

xlabel( 'Time (s) ',' FontSize' , 16) ; ylabel( 'Phase Angle, \theta 
(deg)' , 'FontSize' , 20) ; 

legend( '\theta_v_{_s_o_u_r_c_e}' , '\theta_i_{_s_o_u_r_c_e}' ); 

figure( 'name' ,' Phase Angles from t=[0,0.35]sec' ) 
subplot(2,1,1); grid; 
hold on; 

plot (time, v_phase_angle’^180/pi, ' b' , ' linewidth' , 2) ; 
plot (time, i_phase_angle’^180/pi, ' g' , ' linewidth' , 2) ; 
hold off; 

title( 'Resistive Load (85.7 \Omega) ',' FontSize' ,20); 
xlim([0 0.1]); 

xlabel( 'Time (s) ',' FontSize' ,20); ylabel( 'Phase Angle, \theta 
(deg)' , 'FontSize' , 20) ; 

legend( '\theta_v_{_s_o_u_r_c_e}' , '\theta_i_{_s_o_u_r_c_e}' , 'location' , ' 

southeast' ); 

subplot(2,1,2); grid; 

hold on; 

plot (time, v_phase_angle’^180/pi, ' b' , ' linewidth' , 2) ; 
plot (time, i_phase_angle’^180/pi, ' g' , ' linewidth' , 2) ; 
hold off; 

title (' Inductive Load (0.246 H) Added to the Circuit at t = 0.25 
sec' , 'FontSize' ,20); 
xlim([0.25 0.35]); 

xlabel( 'Time (s) ',' FontSize' ,16); ylabel( 'Phase Angle, \theta 
(deg)' , 'FontSize' , 20) ; 

legend( '\theta_v_{_s_o_u_r_c_e}' , '\theta_i_{_s_o_u_r_c_e}' ); 
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figure( 'name'Calculated PF Angle Between Source Voltage and Source 
Current' ) 

plot (time^ calculated_PF_angle’^180/pi^ ' b' , ' linewidth' , 2 ) } 
xlabel( 'Time (s)' , 'FontSize' ,20); ylabel( 'Calculated PF Angle, \phi 
(deg)' , 'FontSize' , 20) ; 
grid on; xlim([0 tstop]); 

figure( 'name' ,' Corrected and Filtered Power Factor Angle') 

subplot (2,1, 1) , plot (time, corrected_PF_angle’^180/pi, ' b' , ' linewidth' , 2) ; 

grid; xlim([0 tstop]); 

xlabel( 'Time (s)' , 'FontSize' ,20); ylabel( 'Corrected PF Angle, \phi 
(deg)' , 'FontSize' ,20); 

%legend('Unfiltered \phi (\theta_v_{_s_o_u_r_c_e} - 
\theta_i_{_s_o_u_r_c_e})','location','south'); 

subplot (2,1,2) , plot (time, f litered_PF_angle’^180/pi, ' b' , ' linewidth' , 2) ; 
grid; xlim([0 tstop]); 

xlabel( 'Time (s) ',' FontSize' ,20); ylabel( 'Filtered PF Angle, \phi 
(deg)' , 'FontSize' ,20); 

%legend('Filtered \phi (\theta_v_{_s_o_u_r_c_e} - 
\theta_i_{_s_o_u_r_c_e})','location','south'); 
grid on; 

figure( 'name' ,' Power Factor Angle and EMS Current Amplitude') 
subplot (2,1, 1) , plot (time, f iltered_PF_angle’^180/pi, ' b' , ' linewidth' , 2) ; 
grid; xlim([0 tstop]); 

xlabel( 'Time (s) ',' FontSize' ,20); ylabel( 'Filtered PF Angle, \phi 
(deg)' , 'FontSize' ,20); 

%legend('Filtered \phi (\theta_v_{_s_o_u_r_c_e} - 
\theta_i_{_s_o_u_r_c_e})','location','south'); 

subplot(2,1,2),plot(time,i_ems_amplitude, 'g' , 'linewidth' ,2); grid; 
xlim([0 tstop]); 

xlabel( 'Time (s )',' FontSize' , 20) ; ylabel( 'EMS Current Amplitude, 
I^^_e_m_s (A)' , 'FontSize' ,20); 

grid on; 

C. EXPERIMENT OUTPUT PLOT FILE 

% Michael Prato 

% Reactive Power Compensation Using an EMS 
% Oscilloscope Output Plot File 

Clc; 

clear all; 
close all; 

data_vac=xlsread( 'Tek_CHl_Wfm.csv' ); 
data_iems=xlsread( 'Tek_CH2_Wfm.csv' ); 
data_isource=xlsread( 'Tek_CH3_Wfm.csv' ); 
data_iload =xlsread( 'Tek_CH4_Wfm.csv' ); 
len = length(data_vac); 

time_vec=data_vac(15:len, 1) ; 
vac=data_vac(15:len, 2 ) ; 
iems=data_iems(15:len, 2 ) ; 
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iload=data_iload(15:len,2); 
isource=data_isource(15:len^ 2); 
tmin=time_vec(1); 
tmax=max(time_vec); 

data_vac2=xlsread( 'Tek_CHl_Wfmb.csv' ); 
data_iems2=xlsread( 'Tek_CH2_Wfmb.csv' ); 
data_isource2=xlsread( 'Tek_CH3_Wfmb.csv' ); 
data_iload2 =xlsread( 'Tek_CH4_Wfmb.csv' ); 
len = length(data_vac2); 

time_vec2=data_vac2(15:len^1) ; 
vac2=data_vac2(15:len, 2 ) ; 
iems2=data_iems2(15:len, 2) ; 
iload2=data_iload2(15:len, 2) ; 
isource2=data_isource2(15:len, 2) ; 
tmin2=time_vec(1); 
tmax2=max(time_vec) ; 

figure(1); 

subplot(2,1,1), plot(time_vec,vac/50); 
hold on; 

plot(time_vec,isource, ' r' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e 

(A) ' , 'Location' , 'Northeast' ) ; grid on; 
axis([tmin tmax -4 4]); 

subplot(2,1,2),plot(time_vec,iems, 'g' ,time_vec,iload, 'm' ); 
xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 
grid on; 

legend( '-i_e_m_s (A)' , 'i_l_o_a_d (A) ',' Location' ,' Northeast ') ; 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ) ; 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figurel' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figurel' ) ; 

figure(2); 

subplot(2,1,1), plot(time_vec2,vac2/50); 
hold on; 

plot(time_vec2,isource2, ' r' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e 

(A)' , 'Location' , 'Northeast' ) ; 
grid on; 

axis([tmin tmax -4 4]); 
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subplot { 2 , 2 ) ^ plot(time_vec2,iems2 ^ 'g ' a time_vec2,iload2, 'm' ) 
xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 
grid on; 

legend( '-i_e_m_s (A)' ^ 'i_l_o_a_d (A)' ^ 'Location' ^ 'Northeast' ) 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ) ; 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figure2' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figure2' ) ; 

figure(3); 
hold on; 

plot(time_vec2,vac/50) ; 
plot(time_vec2,isource, ' r' ) ; 
plot(time_vec2,iems, ' g' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '- i_e_m_s 

(A)' , 'Location' , 'Northeast' ); 
grid on; 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ) ; 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figure3' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figure3' ) ; 

figure(4); 
hold on; 

plot(time_vec2,vac2/50) ; 
plot(time_vec2,isource2, ' r' ) ; 
plot(time_vec2,iems2, ' g' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '- i_e_m_s 

(A)' , 'Location' , 'Northeast' ) ; 
grid on; 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ) ; 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figure4' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figure4' ) ; 

figure(5); 
hold on; 






plot(time_vec,vac/50) ; 
plot(time_vec,isource, ' r' ) ; 
plot(time_vec,iems^ ' g' ) ; 
plot(time_vec,iload, 'm' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' ^ 'i_s_o_u_r_c_e (A)' , '-i_e_m_s 

(A) ' , ' i_l_o_a_d (A) ' ' Location' , ' Northeast' ) ; 

grid on; 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ) ; 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figure5' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figure3' ) ; 

figure(6); 
hold on; 

plot(time_vec,vac2/50) ; 
plot(time_vec,isource2, ' r' ) ; 
plot(time_vec,iems2, ' g' ) ; 
plot(time_vec,iload2, 'm' ) ; 
hold off; 

xlabel( 'Time (s)' ); 
ylabel( 'Amplitude' ) ; 

legend( 'v_s_o_u_r_c_e/50 (V)' , 'i_s_o_u_r_c_e (A)' , '-i_e_m_s 

(A)' , 'i_l_o_a_d (A)' , 'Location' , 'Northeast' ); 
grid on; 

axis([tmin tmax -4 4]); 

templ=get(gcf, 'position' ); 

temp2=templ+[-500, 0, tempi (3)’^l. 5, 0 ]; 

set(gcf, 'position' ,temp2); 

saveas(gcf, 'figure6' , 'tif' ); 

print(gcf, '-dtiff' , '-r350' , 'figure4' ) ; 
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